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ABSTRACT 

We report the results of XMM-Newton observations of the X-ray bright, radio-loud blazar 
GB B1428H-4217 at a redshift of z = 4.72. We confirm the presence of soft X-ray spectral 
flattening at energies < 0.7 kcV as reported in previous ROSAT and BeppoSAX observations. 
At hard X-ray energies the spectrum is consistent with a power-law although we find the spec- 
tral slope varied between both XMM-Newton observations and is also significantly different 
from that reported previously. Whilst we cannot rule-out intrinsic cold absorption to explain 
the spectral depression, we favour a dust-free warm absorber. Cold absorption requires a col- 
umn density ~ 1.4 — 1.6 x 10^^ cm^^ whilst a warm absorber could have up to ~ 10^'^ cm~^ 
and an ionization parameter ^ 10^. The spectrum of GB B1428H-4217 shows remarkable par- 
allels with that of the z = 4.4 blazar PMN J0525-3343, in which the available evidence is 
also most consistent with a warm absorber model. 
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1 INTRODUCTION 

High-redshift quasars are important objects in the study of the early 
evolution of massive black holes and the interaction with their host 
galaxies. A number of radio-loud quasars have been found to show 
' the spectral characteristics typical of blazars . GB B1428+4217 
! at z = 4.72 (Hook & McMahon 1998; Fabia n et alJll997lll998l 
11999') displays spectral flattening at soft X-ray energies in ROSAT 
(Boiler et al. 2000) and BeppoSAX {Fahian et al. 2001a) data. Here 
we confirm the flattening with XMM-Newton observations. Soft 
' X-ray flattening h as also been reported f rom RX J1028. 6-08 44 
' jYuan et aljEoOO) and PMN 10525-3343 jpabian et alj|200ld) at 
redshifts of z = 4.28 and z — 4.4 respectively. XMM-Newton 
data shows only marginal evid ence f or flattening in the case of 
RX J1028.6-0844 iGruDe et al] |2004. Conversely, the ef fect has 
been confirmed in PMN J0525-3343 <Worslev et aT]|2004ll) . Spec- 
tral flattening is consistent with a trend seen in a number of stud- 
ies of low er-redshift radio-loud quasars (e.g. ICappj etaljfl997t 
iFiore et al.. 1998:.Reeves & Turner 2000) . 

The most consistent explanation for the flattening is intrinsic 
absorption with column densities of the order iq'^'^~'^^ cm~^. In- 
tergalactic absorption systems remain an unlikely possibility given 
their low metallicities. Furthermore, in order to account for the flat- 
tening the line of sight to the blazar would need to intercept two or 
more very high column density systems of which the probability is 
very low. Certainly this explanation cannot account for the spectral 
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depressions seen in an increasing number of objects. An alternative 
explanation for flattening is a break in the underlying blazar con- 
tinuum. Hard X-ray production in blazars is likely to be due to the 
inverse Compton scattering of seed photons by relativistic electrons 
in the jet. A spectral break can arise through a low energy cut-off 
in the electron population or by a peaked seed photon spectrum, 
however the break appears much too sharp for this to be plausible. 

Given an intrinsic absorber hypothesis there remains much 
discussion about the nature of the absorber. X-ray data alone have 
been unsuccessful in breaking the degeneracy between a cold (neu- 
tral) absorber or one that is warm (ionized). Optical observations 
however, can play a crucial role in this regard. A number of the 
blazars showing soft X-ray flattening have optical spectra (corre- 
sponding to rest frame UV emission) which indicate very little dust 
along the line of sight or show any evidence of a Lyman-limit sys- 
tem (for the optical data f or GB B 14 28+4217 a nd PMN J0525- 
3343 see lHook & McMahon 1998 and lPerouxet aL.200L) . These 
point to a dust-free yet metal-enriched absorber where hydro- 
gen has been ionized, consistent with the warm absorber such as 
that s een in Seyfert galaxies (e.g. Reynolds 1997; Crenshaw et al] 
I1999I) . Strong intrinsic C IV absorption is common in these galaxies 
and it is also a feature of PMN J0525-3343. 



2 PREVIOUS OBSERVATIONS 

GB B 1428-1-4217 was discovered bv lHook & McMahoij <1998h in 
an optical sur vey for high-redshift quasars. An X-ray follow-up was 
conducted bv lFabian etai] ll99H) with the ROSAT HRI and A5CA 
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Table 1. Details of XMM-Newton observations. 



Revolution 


Date 


Good Exposure Times (s) 






PN MOS-1 MOS-2 


549 


2002 Dec 09 


2625 4055 4260 


569 


2003 Jan 17 


11534 14241 14235 



jpabian et alll998h . Spectral fits were consistent with a power-law 
continuum with a photon index F = 1.29 although the ASCA SIS 
data alone were suggestive of some excess absorption. Subsequent 
ROSAT observations showed that X-ray flux varied by a factor of 
two over a time-scale of around 2.5 days in the rest frame of the 
source. Radio data also shows variability, with amplitudes of ~ 
40/15 per cent on rest frame time-scales of weeks/days. 

More evidence of X-ray variability was seen in a/?05ArPSPC 
observation ( Boiler et al. 200C), where changes of ~ 25 per cent 
occur on rest frame timescales of 1.8 hours. The low-energy sensi- 
tivity of the PSPC detector clearly revealed a spectral depression in 
the soft X-ray band, consistent with an absorbing column density 
of A'h ~ 1.5 X 10^^ cm~^ (if intrinsic to the source). 

Spectral flatten ing was further confirmed by BeppoSAX 
jpabian et aTlbOOlalh . The fitted spectrum was that of a T = 1.45 
power-law with an intrinsic absorption column density of A'h ~ 
8 X 10^^ cm~^. Warm absorption models were considered, with 
column densities of order A'^h ~ 10^^ '^ cm~^ and ionization pa- 
rameters of order ^ ^ 10^ being marginally consistent with results 
from BeppoSAX as well as those from the previous ROSAT and 
ASCA observations. 



3 XMM-A'£iyroA' OBSERVATIONS 

3.1 Data reduction 

XMM-Newton observed GB B 1428-1-4217 on 2002 December 9 and 
again on 2003 January 17 (Table 0. All observations were per- 
formed in full frame imaging mode with the thin filter. The EPIC 
data reduction was performed with the Scientific Analysis Software 
(SAS) V5.4.1. Periods of background flaring were excluded and the 
data were filtered in the standard way to include only the recom- 
mended events for spectral analysis (Ehle et al. 2003) i.e. using the 
event pattern ranges — 4 and — 12 plus the 'FLAG=0' and 
'#XMMEA^M' filters for the PN and MOS respectively. 

Spectra were extracted from circular regions centred on the 
source of approximate radius 40" and 60" for the PN and MOS re- 
spectively. Backgrounds were ex tracted from nearb y source free re- 
gions following the guidelines in lEhle et a l. ( 2003"). SAS tasks AR- 
FGEN and RMFGEN were used to generate the appropriate response 
files. Background subtracted spectra were produced and grouped 
to a minimum of 20 counts in each bin. Spectral analysis was per- 
formed using XSPEC Vl 1.3.0 lArnaud, 1996.1 . 

3.2 Spectral properties 

Spectral fits to a simple power-law in clusive of neutral galactic ab- 
sorption (A^H = 1-40 X 10^" cm~^: lEivis et al.i.1994) were made 
over the energy range 1 — 12 keV. The PN and MOS data were 
fitted simultaneously, but the revolution 549 and 569 observations 
were treated separately due to source variability. 1 — 12 keV power- 
law photon indices of T = 1.88 ± 0.05 and V = 1.73 ± 0.03 were 
obtained, with 2 - lOkeV fluxes of 1.22 x 10"^^ erg cm"^ 
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Figure 1. The upper panel shows the ratio of the revolution 549 data to 
a simple F = 1.88 power-law (including Galactic absorption of A'h = 
1.40 X 10^" cm~^; fitted over 1 — 12keV). The lower panel shows the 
revolution 569 data to a F = 1.73 power-law (with the same Galactic cor- 
rection and fit range). The PN data are shown in both instances as black 
crosses with the MOS-l/MOS-2 data as black/grey diamonds respectively. 

and 1.09 x lO"^'^ erg cm~'^ respectively. Reduced chi-squared 
values for the fits were = 0.819 and = 0.883 (with v = 100 
and u = 283 degrees of freedom) respectively. These can be com- 
pared to only = F06 (u = 385) with an attempt to fit the same 
model to the data from both observations simultaneously. 

A significant depression in the soft X-ray spectrum is seen 
when the whole 0.3 — 12 keV energy range is considered. The chi- 
squared statistic of the power-law models rises to xt = 1-245 and 
xt ~ 1.302 (for u = 173 and u = 430) respectively. Fig.Qshows 
the ratio of the revolution 549 and 569 data to the 1 — 12 keV fitted 
power-law. Spectral flattening is clear in both observations, with a 
break from the power-law occurring at energies ~ 0.6 — 1 keV and 
the data/model ratio decreasing to ~ 0.7 at < 0.4 keV. 

Even though there has been a significant decrease in flux be- 
tween the observations (and possibly spectral slope variations) it 
remains instructive to consider the combined PN spectra in the 
search for spectral features which may exist at the same energy 
independently of other variations. To this end the PN camera good 
events from both observations were merged and a single spectrum 
extracted. Fig. |2| shows the ratio of the data to a 1 — 12 keV fitted 
power-law with F — 1.77 (plotted in the rest energy frame of the 
blazar). Apart from the distinctive spectral flattening there appear 
to be no statistically significant features. 

The spectral break can be modelled in a number of ways. Ta- 
ble |2| shows the results of broken power-law and power-law plus 
absorption models. The best-fitting local absorption column den- 
sity exceeds the expected galactic level at the 7-sigma level. If the 
excess is modelled as cold absorption at the source then column 
densities of 1.4 — 1.6 x lO'^^ cm^^ are required. Fig.|3|shows the 
confidence contours in the photon index and column density (for 
intrinsic cold absorption) for the two observations. 

3.3 Variability 

Background-corrected lightcurves were extracted for both obser- 
vations. The good- time-interval lightcurve from the revolution 549 
observation is fully consistent with a constant flux (the PN count- 
rate is 0.73 cts corresponding to 1.22 x 10~^^ erg cm~^ s~^). 
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Table 2. Summary of the spectral fits performed to the revolution 549 and 549 data. Fits were made in each case, to the PN and MOS spectra jointly over the 
energy range 0.3 — 12 keV. Errors indicated are 1-sigma for one interesting parameter. 



Model 



local A^H intrinsic TVjj 

(1020cm-2) (1022 cm~2) 



-^brcak 

(keV) 



Rev. 549 data: 

Broken power-law 1.40 (Galactic) 
power-law, free local absorption 5.0 ±0.5 

power-law, free intrinsic absorption 1.40 (Galactic) 1.4 it 0.4 



1.1 ±0.3 0.7 ±0.2 



1.83 ±0.03 
1.92 ±0.05 
1.87 ±0.04 



0.851 171 
0.840 172 
0.844 172 



Rev. 569 data: 

Broken power-law 1.40 (Galactic) 
power-law, free local absorption 5.1 ±0.5 

power-law, free intrinsic absorption 1.40 (Galactic) 1.6 ±0.2 



0.7±0.2 0.6±0.1 



1.71 ±0.02 
1.79 ±0.03 
1.76 ±0.02 



0.932 428 
0.949 429 
0.931 429 




2 5 10 20 50 

Rest frame energy (keV) 

Figure 2. The combined PN camera data of both observations, plotted as 
a ratio to a fitted F = 1.77 power-law (including Galactic absorption of 
A^H = 1-40 X 10^" cm~2) xhe energy has been adjusted to the source 
rest frame; the observed range is identical to that covered in Fig.lTI 




Figure 3. Confidence contours (68, 90 and 99 per cent) in the photon index 
and absorber column density for a power-law plus intrinsic cold absorber 
model (also including galactic absorption). The separate contours for revo- 
lutions 549 and 569 are shown in grey and black respectively. 



Unfortunately the remaining ~ 60 per cent of the observation 
was affected by solar flare activity and rehable background sub- 
traction could not be performed. The revolution 549 observation 
suffered from no flare periods and ~ 12 ks of PN data were ex- 
tracted, again this is consistent (at the 90 per cent level) with a 
constant flux (a PN count-rate of 0.54 ctss"^ corresponding to 
1.09x lO"'^'^ erg cm~^ s~^). The ~ 10 per cent reduction between 
the observations corresponds to a variability timescale of ~ 7 days 
in the rest frame. 

In addition to the decrease in luminosity there are differences 
in the spectral properties of the source between observations. The 
revolution 549 data are well fit (in the 1 — 12keV range) by a 
simple power-law with a photon index of F ~ 1.9 whilst this is 
hardened to F < 1.8 during the revolution 569 observation. The 
break energy and soft X-ray photon index (and equivalently, the 
column densities in the absorption models) are in good agreement 
between the observations. 

There would also seem to be significant variations in the spec- 
tral shape between the XMM-Newton and earlier observations. The 
continuum slopes (i.e. after correction for any soft X-ray flattening) 
reporte d by ROSAT and BeppoSAX l Boiler et al. 2000 ; Fabian et al. 
^ooTP) were 1.4 ± 0.2 and 1.45 ± 0.1 respectively, considerably 
harder than the ~ 1.7— 1.9 found with XMM-Newton. The spectral 
break is characterised by the column density required of an intrin- 
sic cold absorber, the ROSAT and BeppoSAX values of which are 



1.5 ± 0.3 X 10^^ cm"^ and 7.8t.l'l x 10^^ cm"^, consistent with 
each other and the ~ 1.4 — 1.6 x 10^^ cm~^ found here. 

There is marginal evidence of curvature in the residuals of the 
XMM-Newton fits to the 1 — 12 keV power-law, with excesses at 
~ 1 — 2 keV and ~ 5 — 10 keV and a deficit over ~ 2 — 5 keV 
(Fig.0. The > 5keV range appears to have a harder power-law 
(F ~ 1.5), which is in better agreement with the ROSAT and Bep- 
poSAX results. A steeper 1 — 12 keV slope and curved residuals 
in the XMM-Newton data could possibly then be the result of a 
~ 1 — 3 keV variable excess (~ 5 — 20 keV rest frame) super- 
imposed upon the soft X-ray break. 

To investigate further a difference spectrum was formed by 
subtracting the revolution 569 spectrum from the revolution 549 
one. The difference is well-fitted by a power-law model (inclu- 
sive of Galactic absorption and intrinsic cold absorption of A'h = 
1.5 X 10^^ cm"^) with a spectral slope of T = 2.3 ± 0.3. A two 
power-law model (again inclusive of Galactic and intrinsic absorp- 
tion) is also an acceptable fit to both data sets. The soft power-law 
component has F ~ 1.8 — 2.6 and the hard power-law has F < 1.7. 

The curved residuals, 1 — 3keV excess and the success of 
the two power-law model point to the intriguing possibility that 
the spectrum may contain a additional component. The steepness 
(F ~ 2.3) of the soft component is consistent with the high-energy 
cut-off of a second inverse-Compton contribution. Possibilities for 
such a component are an underlying pile-up in the EC seed photon 
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distribution (see l4.2> . an EC component witii a different external 
photon source, SSC emission or a bulk Compton component due to 
inverse-Compton scattering from a cold electron population. 



4 DISCUSSION 

There are two main possibilities for the origin of the spectral flat- 
tening observed in GB B 1428+4217. The first, and our favoured 
explanation, is that the flattening is due to excess intrinsic absorp- 
tion. The second hypothesis attributes the effect to a break in the 
underlying blazar emission continuum. These two alternatives are 
now discussed separately. 

4.1 Absorption 

A galactic origin to the excess column density is extremely un- 
likely; it would require a factor ~ 3.6 increase in the galactic col- 
umn density, highly localised in the direction of GB B 1428-1-4217. 
Intergalactic absorption is also highly unlikely, the only possibil- 
ity being a line of sight that passes through two or more very high 
column d ensity damped Ly-g systems for which the probab ility is 
very low I'O'Flahertv & Jakobserj |l997tlZwaan et m h99'^. If an 
absorber is the origin of the X-ray flattening it must therefore be 
metal-rich and intrinsic to the source. A cold absorber model re- 
quires a column density in the region 1.4 — 1.6 x ICp'^ cm~^. 

A warm (ionized) absorber is able to reproduce both low 
UV and high X-ray column densities. Photoionized gas could pro- 
duce significant photoelectric absorption in the soft X-ray band 
(mainly due to oxygen and neon species) whilst remaining transpar- 
ent at UV wavelengths (hydrogen/helium both ionized). Addition- 
ally, any dust grains could be destroyed by the intense photon-flux, 
resulting in low levels of extinction in the UV band. 

The photoionization code CLOUDY ( Ferla nd et alJll99^ was 
used to model the warm absorber hypothesis. Fig. |4] shows the 
confidence contours in the warm absorber ionization parameter 
and column density for the two separate observations. Neither set 
of contours strongly constrains the ionization parameter provided 
that ^ < 10^ erg cm s~^. The column density is in the range 
10^^ < A^H < 10^^ cm-2. 

It is difficult to break this degeneracy with X-ray observa- 
tions: The conclusive signature of a warm absorber would be the 
detection of photoelectric absorption edges, principally those aris- 
ing from 0/Ne. Due to the high-z of the source however, these fea- 
tures are redshifted out of the observed band to energies < 0.2 keV. 
A strong discriminator would be evidence of Lyman-limit absorp- 
tion at 5217 A (the present optical data start at 5970 A). This would 
place constraints on the neutral hydrogen column density and thus 
on the ionization state of the absorber. The lack of an optically- 
thick Lyman-limit system would immediately place an upper limit 
to the H I column density of ~ 10^^ cm^^, requiring the warm ab- 
sorption model to have an ionization parameter 5 > 1 erg cm s^^. 

4.2 Intrinsic spectral properties 

The revolution 549 and 569 spectra can be acceptably fitted with a 
broken power-law model with a break at an energy ~ 0.6 — 0.7 keV 
down to a flat spectral index ~ 0.7 — 1.1. In the rest frame (Fig.|2j 
the break energy is ~ 4 keV and is very distinct, occurring over a 
range of ~ 2 keV. The z = 4.4 blazar PMN J0525~3343 is similar 
with a sharp spectral break at a rest frame energy of ~ 4.5 keV 
over an energy range < 5 keV. 




N-3 _g _i 1 2 3 



Log f (erg cm s ^) 

Figure 4. Confidence contours (68, 90 and 99 per cent) in the photon index 
and column density for a power-law plus intrinsic cold absorber model (also 
including galactic absorption). The separate contours for the revolution 549 
and 569 observations are shown in grey and black respectively. 

The dominant emission components in blazars are beamed 
synchrotron radiation from the jet (peaking in the infra-red to soft 
X-ray band) and a high-energy component (dominating the hard 
X-ray to 7-ray regimes). The high-energy emission arises through 
the inverse Comptonization of soft photons by highly relativistic 
electrons in the jet plasma. Sources of the soft photons are the nu- 
clear opticalAJV emission or the local synchrotron radiation within 
the jet. For the brightest blazars (such as GB B1428-I-4217 and 
PMN J0525-3343) the fi rst, known as the 'e xternal Compton' (EC) 
mechanism, is dominant JSikora et alll994l) . 

A flattening in the spectrum could be expected if there was a 
low energy cut-off in the energy distribution of the electron pop- 
ulation in the jet. In the case of the EC process it is also possible 
to produce a spectral break by requiring the seed photon distribu- 
tion to be sharply peaked at and/or rapidly decline below a cer- 
tain frequency. There are several problems with these explanations: 
Firstly, there is no evidence of any spectral breaks in the spectra 
of nearby radio-loud quasars, and secondly, there seems to be no 
plausible way to account for the sharpness of the break (only a few 
keV in the rest frame). For a m ore detailed d i scussio n of the break- 
producing mechanisms see e.g. lFabianetalJ<200ld) . 

4.3 Comparisons with other objects 

Fig. |5| shows the rest frame combined PN spectrum of 
GB B 1428-1-4217 (as a ratio to the power-law above the spec- 
tral break) along with the equivalent result for PMN J0525-3343 
(Worslev et al. 2004). The agreement in the position and structure 
of the break is remarkable. Since the plots are of ratios to the model 
the instrumental response (and Galactic absorption) have been re- 
moved. In both cases the break from the power-law occurs over a 
range of < 5 keV in the source frame and to a photon index of 
r ~ 1. Both also show a steeper drop in the flux at 3 keV. 

If the spectral flattening is indeed due to warm absorption then 
the absorbers for GB BI428-H42I7 and PMN J0525-3343 are al- 
most identical in column density (of ~ 2 — 4 x 10^^ cm~'^). There 
is also evidence for a similar column in the XMM-Newton observa- 
tions of RX J1028.6-0844 I Grupe et al. 2004) where joint fits to the 
XMM-Newton and previous ASCA data iYuan et all2000l) are con- 
sistent with an excess intrinsic column density of ~ 1 x 10^^ cm~^ . 
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Figure 5. The ratio of combined PN camera spectra to a liard-band fit- 
ted power-law. The filled squares show the data for GB B 1428-1-4217; the 
power-law has F = 1.77 and was fitted over the observed 1 — 12keV 
range. The open circles show the data for PMN J0525-3343; the power-law 
has r = 1.71 and was fitted over the observed 2 — 12keV range. The 
spectra are plotted in their respective rest frame energies. The spectra have 
also been renormalised to a data/model ratio of unity over the 5 — 10 keV 
rest frame range. 

5 CONCLUSIONS 

We have studied the X-ray spectnam of the high-redshift, radio-loud 
blazar GB B 1428-1-4217, providing definitive confirmation of the 
soft X-ray spectral flattening as previously observed with ROSAT 
and BeppoSAX. The flattening is not consistent with Galactic or 
intergalactic absorption and we propose that the effect is due to an 
absorber that is intrinsic to the source, or to an underlying break in 
the blazar continuum. 

A continuum break is highly speculative given the present 
lack of knowledge concerning emission processes in the nuclear re- 
gions. Continuum breaks are not seen in nearby radio-loud objects 
where they should be readily detectable and there is no satisfactory 
way to account for the sharpness of the break seen. An intrinsic 
absorber is therefore our favoured explanation. Cold absorption re- 
quires a column density A'h ^ 1-4 — 1.6 x lO^'^ cm~^, whilst a 
warm absorber could have a column density of up to ~ 10^^ cm~'^ 
and an ionization parameter ^ ~ 10^. X-ray data alone are unable 
to distinguish between a warm or cold absorber model. Good con- 
straints are expected from future optical observations which could 
place strong limits on the neutral hydrogen column density. 

It is interesti ng to note the p arallels with the blazar 
PMN J0525-3343 iWorslev et alJl2004^ i.e. a relatively blue UV 
continuum with strong Ly-a and C IV emission lines as well as the 
flat hard X-ray spectral slope and the remarkably similar form of 
the soft X-ray depression. The PMN J0525-3343 X-ray data are 
consistent with both warm/cold absorption but optical observations 
rule out a cold absorber scenario. The lack of significant redden- 
ing in the UV spectra of both objects is consistent with a dust-free 
warm absorber. 

Comparisons can also be drawn to the z = 3.27 radio- 
loud quasar PKS 2126-0158 which also shows definite so ft X-ray 
absorption <Ferrero & Brinkmand Eoosl iFiore et alj 12003 ). Both 
waim and cold absorbers are consistent although ionized absorp- 
tion is again most able to explain the lack of significant UV red- 
dening. The implied iron abundance of a warm absorber however 
is at odds with the super-solar abundances implied by broad emis- 
sion lines. 



The spectrum of GB B 1428-1-4217 are suggestive of an ad- 
ditional component contributing to an excess of emission in the 
~ 5 — 10 keV rest-frame energy range. There are a number of pos- 
sibilities for the emission including a pile-up in the electron energy 
distribution or a feature caused by bulk Comptonization, however 
the detection is marginal and further work is required. 
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